Chapter  1 

Network  Abstract  Linear  Programming  with 
Application  to  Cooperative  Target  Localization* * 


Giuseppe  Notarstefano  and  Francesco  Bullo 


Abstract  We  identify  a  novel  class  of  distributed  optimization  problems,  namely  a 
networked  version  of  abstract  linear  programming.  For  such  problems  we  propose 
distributed  algorithms  for  networks  with  various  connectivity  and/or  memory  con¬ 
straints.  Finally,  we  show  how  a  suitable  target  localization  problem  can  be  tackled 
through  appropriate  linear  programs. 


1.1  Introduction 

This  paper  focuses  on  a  class  of  distributed  computing  problems  and  on  its  appli¬ 
cations  to  cooperative  target  localization  in  sensor  networks.  To  do  so,  we  study 
abstract  linear  programming,  that  is,  a  generalized  version  of  linear  programming 
that  was  introduced  by  Matousek,  Sharir  and  Welzl  in  [1]  and  extended  by  Gartner 
in  [2].  Abstract  linear  programming  is  applicable  also  to  some  geometric  optimiza¬ 
tion  problems,  such  as  the  minimum  enclosing  ball,  the  minimum  enclosing  stripe 
and  the  minimum  enclosing  annulus.  These  geometric  optimization  problems  are 
relevant  in  the  design  of  efficient  robotic  algorithms  for  minimum-time  formation 
control  problems  as  shown  in  [3]. 
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Linear  programming  and  its  generalizations  have  received  widespread  attention 
in  the  literature.  The  following  references  are  most  relevant  in  our  treatment.  The 
earliest  (deterministic)  algorithm  that  solves  a  linear  program  in  a  fixed  number  of 
variables  subject  to  n  linear  inequalities  in  time  O(n)  is  given  in  [4].  An  efficient  ran¬ 
domized  incremental  algorithm  for  linear  programming  is  proposed  in  [1],  where  a 
linear  program  in  d  variables  subject  to  n  linear  inequalities  is  solved  in  expected 
time  0{d2n  +  e°^‘no^)\  the  expectation  is  taken  over  the  internal  randomizations 
executed  by  the  algorithm.  An  elegant  survey  on  randomized  methods  in  linear  pro¬ 
gramming  is  [5];  see  also  [6].  The  survey  [7],  see  also  [8],  discusses  the  application 
of  abstract  linear  programming  to  a  number  of  geometric  optimization  problems. 
Regarding  parallel  computation  approaches  to  linear  programming,  we  only  note 
that  linear  programs  with  n  linear  inequalities  can  be  solved  [9]  by  n  parallel  pro¬ 
cessors  in  time  0(  (log  log  (w)  )rf).  The  approach  in  [9]  and  the  ones  in  the  references 
therein  are,  however,  limited  to  parallel  random-access  machines  (usually  denoted 
PRAM),  where  a  shared  memory  is  readable  and  writable  to  all  processors.  In  this 
paper,  we  focus  on  networks  described  by  arbitrary  graphs. 

The  problem  of  target  localization  has  been  widely  investigated  and  the  related 
literature  is  therefore  quite  rich.  Recently,  the  interest  in  sensor  networks  and  dis¬ 
tributed  computation  has  lightened  the  attention  on  this  problem  from  this  new  per¬ 
spective.  A  good  reference  for  localization  and  tracking  in  sensor  networks  is  [10]. 
Two  approaches  may  be  used  to  tackle  target  localization,  a  stochastic  and  a  deter¬ 
ministic  one.  As  regards  the  deterministic  approach  (the  one  we  use  in  this  paper), 
a  set  membership  estimation  technique  was  proposed  in  [11].  Recently,  a  sensor 
selection  problem  for  target  tracking  was  studied  in  [12].  References  for  target  lo¬ 
calization  and  tracking  in  sensor  networks,  even  by  use  of  a  stochastic  approach, 
may  be  found  therein. 

The  contributions  of  this  paper  are  three-fold.  First,  we  identify  a  class  of  dis¬ 
tributed  optimization  problems  that  appears  to  be  novel  and  of  intrinsic  interest. 
Second,  we  propose  a  novel  simple  algorithmic  methodology  to  solve  these  prob¬ 
lems  in  networks  with  various  connectivity  and/or  memory  constraints.  Specifically, 
we  propose  three  algorithms,  prove  their  correctness  and  establish  halting  condi¬ 
tions.  Finally,  we  illustrate  how  these  distributed  computation  problems  are  relevant 
for  distributed  target  localization  in  sensor  networks.  Specifically,  we  cast  the  tar¬ 
get  localization  problem  in  the  problem  of  approximating  the  intersection  of  convex 
poly  topes  by  using  a  small  number  of  halfplanes  (among  all  those  defining  the  poly¬ 
topes).  We  show  that  suitable  linear  programs  running  in  parallel,  in  fact,  solve  the 
problem,  so  that  the  proposed  distributed  algorithms  may  be  used. 

The  paper  is  organized  as  follows.  Section  1.2  introduces  abstract  linear  pro¬ 
grams.  Section  1.3  introduces  network  models.  Section  1.4  contains  the  definition 
of  network  abstract  linear  programs  and  the  proposed  distributed  algorithms.  Sec¬ 
tion  1.5  shows  the  relevance  of  the  proposed  distributed  computing  algorithms  in 
the  context  of  cooperative  target  localization. 
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Notation 


We  let  N,  No,  and  R+  denote  the  natural  numbers,  the  non-negative  integer  numbers, 
and  the  positive  real  numbers,  respectively.  For  f,g  :  N  — >  R,  we  say  that  /  £  0(g) 
if  there  exist  no  £  N  and  k  £  M+  such  that  \f(n)\  <  k\g(n)  \  for  all  n  >  no. 


1.2  Abstract  linear  programming 

In  this  section  we  present  an  abstract  framework  that  captures  a  wide  class  of  opti¬ 
mization  problems  including  linear  programming  and  various  geometric  optimiza¬ 
tion  problems.  These  problems  are  known  as  abstract  linear  programs  (or  LP-type 
problems).  They  can  be  considered  a  generalization  of  linear  programming  in  the 
sense  that  they  share  some  important  properties.  A  comprehensive  analysis  of  these 
problems  may  be  found  for  example  in  [7], 


1.2.1  Abstract  framework 

We  consider  optimization  problems  specified  by  a  pair  (H,  ft)),  where  H  is  a  finite 
set,  and  ft) :  2H  — >  £2  is  a  function  with  values  in  a  linearly  ordered  set  (£2 ,  <);  we 
assume  that  £2  has  a  minimum  value  —  The  elements  of  H  are  called  constraints, 
and  for  G  CH,  Co(G)  is  called  the  value  of  G.  Intuitively,  ft)(G)  is  the  smallest  value 
attainable  by  a  certain  objective  function  while  satisfying  the  constraints  of  G.  An 
optimization  problem  of  this  sort  is  called  abstract  linear  program  if  the  following 
two  axioms  are  satisfied: 

(i)  Monotonicity:  if  F  C  G  C  H,  then  co(F)  <  Co(G); 

(ii)  Locality:  if  F  C  G  C  Ft  with  —  °°  <  (o(F)  =  ft)(G),  then,  for  all  h  £  H, 


A  set  B  C  H  is  minimal  if  co(B)  >  Co(B')  for  all  proper  subsets  B'  of  B.  A  minimal 
set  B  with  —  °°  <  Co(B)  is  a  basis.  Given  G  C  H,  a  basis  of  G  is  a  minimal  subset 
B  C  G,  such  that  —  °°  <  co(B)  =  ft)(G).  A  constraint  h  is  said  to  be  violated  by  G,  if 
®(G)  <  ft)(GU {/)}). 

The  solution  of  an  abstract  linear  program  (H,  ft))  is  a  minimal  set  Bu  C  H  with 
the  property  that  co(Bh)  =  co(H).  The  combinatorial  dimension  8  of  ( H,co )  is  the 
maximum  cardinality  of  any  basis.  Finally,  an  abstract  linear  program  is  called  basis 
regular  if,  for  any  basis  with  card(B)  =  5  and  any  constraint  h  £  Ft,  every  basis 
of  BU  {/)}  has  the  same  cardinality  of  B.  We  now  define  two  important  primitive 
operations  that  are  useful  to  solve  abstract  linear  programs. 

(i)  Violation  test:  given  a  constraint  h  and  a  basis  B,  it  tests  whether  h  is  violated  by 

B:  we  denote  this  primitive  by  Viol (B,h)\ 
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(ii)  Basis  computation:  given  a  constraint  h  and  a  basis  B,  it  computes  a  basis  of 
BU{li}',  we  denote  this  primitive  by  Basis (B,h). 

Remark  1  (Examples  of  abstract  linear  programs).  We  present  three  useful  geomet¬ 
ric  examples;  see  Figure  1.1. 

(i)  Smallest  enclosing  ball:  Given  n  points  in  Rf/,  compute  the  center  and  radius  of 
the  ball  of  smallest  volume  containing  all  the  points.  This  problem  has  combina¬ 
torial  dimension  d  +  1 . 

(ii)  Smallest  enclosing  stripe:  Given  n  points  in  R2  in  generic  positions,  compute  the 
center  and  the  width  of  the  stripe  of  smallest  width  containing  all  the  points.  This 
problem  has  combinatorial  dimension  5. 

(iii)  Smallest  enclosing  annulus:  Given  n  points  in  R2,  compute  the  center  and  the  two 
radiuses  of  the  annulus  of  smallest  area  containing  all  the  points.  This  problem 
has  combinatorial  dimension  4. 


Fig.  1.1  Smallest  enclosing  ball,  stripe  and  annulus 


More  examples  are  discussed  in  [1,  2,  5,  7]. 


□ 


1.2.2  Randomized  sub-exponential  algorithm 

A  randomized  algorithm  for  solving  abstract  linear  programs  has  been  proposed 
in  [1].  Such  algorithm  has  linear  expected  running  time  in  terms  of  the  number  of 
constraints,  whenever  the  combinatorial  dimension  8  is  fixed,  and  subexponential 
in  8.  The  algorithm,  called  SUBEX_lp,  has  a  recursive  structure  and  is  based  on  the 
two  primitives  introduced  above,  i.e.,  the  violation  test  and  the  basis  computation 
primitives.  For  simplicity,  we  assume  here  that  such  primitives  may  be  implemented 
in  constant  time,  independent  of  the  number  of  constraints.  Given  a  set  of  constraints 
G  and  a  candidate  basis  C  C  G,  the  algorithm  is  as  follows. 


Title  Suppressed  Due  to  Excessive  Length 


5 


function  SUBEX.lp (G,C) 

if  G  =  C,  then  return  C 
else 

choose  a  random  h  £  G\C 
B  ■=  SUBEX_lp(G\{/i},C) 
if  Viol(Z?,/t),  i.e.,  h  is  violated 
by  B, 

return 

SUBEX_lp(G,  Basis  (B,/i)) 

else  return  B 
end  if 
end  if 

For  the  abstract  linear  program  (H .  to),  the  routine  is  invoked  with  SUBEX  Lp(H.B), 
given  any  initial  candidate  basis  B. 

In  [1]  the  expected  completion  time  for  the  SUBEX.lp  algorithm  in  conjunction 
with  Clarkson’s  algorithms  was  shown  to  be  in  0(d2n+eol'^dXogd'l)  for  basis  regular 
abstract  linear  programs.  In  [5]  the  result  was  extended  to  problems  that  are  not  basis 
regular. 


1.3  Network  models 

Following  [13],  we  define  a  synchronous  network  system  as  a  “collection  of  com¬ 
puting  elements  located  at  nodes  of  a  directed  network  graph.”  These  computing 
elements  are  sometimes  called  processors. 


1.3.1  Digraphs  and  connectivity 

We  let  I  =  { 1 , . . . , «}  and  let  (.jl  =  (I.E)  denote  a  directed  graph,  where  I  is  the  set  of 
nodes  and  E  C  /  x  I  is  the  set  of  edges.  For  each  node  i  of  <$ ,  the  number  of  edges 
going  out  from  (coming  into)  node  i  is  called  out-degree  ( in-degree )  and  is  denoted 
outdegW  (indeg^).  The  set  of  outgoing  (incoming)  neighbors  of  node  i  are  the 
set  of  nodes  to  (from)  which  there  are  edges  from  (to)  i.  They  are  denoted  .Aj_)(i) 
and  respectively.  A  directed  graph  is  called  strongly  connected  if,  for  every 

pair  of  nodes  (i.j)  £  I  x  I,  there  exists  a  path  of  directed  edges  that  goes  from  i  to 
j.  In  a  strongly  connected  digraph,  the  minimum  number  of  edges  between  node 
i  and  j  is  called  the  distance  from  i  to  j  and  is  denoted  dist (i,j).  The  maximum 
dist(t,y)  taken  over  all  pairs  (i.j)  is  the  diameter  and  is  denoted  diam(Sf).  Finally, 
we  consider  time-dependent  directed  graphs  of  the  form  1 1— >  '.^(t)  =  ( I,E(t )).  The 
time-dependent  directed  graph  Sf  is  jointly  strongly  connected  if,  for  every  t  £  No, 


yj(jfte$(x)  is  strongly  connected. 
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Moreover,  the  time-dependent  directed  graph  Sf  is  uniformly  strongly  connected  if, 
there  exists  S  >  0  s.t.  for  every  t  £  No 

Uf+^(T)  is  strongly  connected. 


1.3.2  Synchronous  networks  and  distributed  algorithms 

Strictly  speaking,  a  synchronous  network  is  a  directed  graph  =  (/,£cmm)  where 
the  set  I  =  {1, . . .  ,n}  is  the  set  of  identifiers  of  the  computing  elements,  and  the 
time-dependent  map  Ecmm  :  No  — >  2/x/  is  the  communication  edge  map  with  the 
following  property:  an  edge  (i,j)  belongs  to  £cmm(i)  if  and  only  if  processor  i  can 
communicate  to  processor  j  at  time  t. 

Definition  1  (Distributed  algorithm).  Let  &  =  (/,  £cmm)  be  a  synchronous  net¬ 
work.  A  distributed  algorithm  consists  of  the  sets 

-  W,  set  of  “logical”  states  for  all  i  £  /; 

-  Wo  C  W,  subset  of  allowable  initial  values; 

-  M,  message  alphabet,  including  the  null  symbol; 

and  the  maps 

-  msg  :  W  x  I  — >  M,  message-generation  function; 

-  stf :  W  x  M"  — >  W,  state-transition  function.  □ 

Execution  of  the  network  begins  with  all  processors  in  their  start  states  and  all 
channels  empty.  Then  the  processors  repeatedly  perform  the  following  two  actions. 
First,  the  ;th  processor  sends  to  each  of  its  outgoing  neighbors  in  the  communication 
graph  a  message  (possibly  the  null  message)  computed  by  applying  the  message- 
generation  function  to  the  current  value  of  wM.  After  a  negligible  period  of  time, 
the  ;th  processor  computes  the  new  value  of  its  logical  variables  vv-'l  by  applying  the 
state-transition  function  to  the  current  value  of  w^,  and  to  the  incoming  messages 
(present  in  each  communication  edge).  The  combination  of  the  two  actions  is  called 
a  communication  round  or  simply  a  round. 

In  this  execution  scheme  we  have  assumed  that  each  processor  executes  all  the 
calculations  in  one  round.  If  it  is  not  possible  to  upper  bound  the  execution-time  of 
the  algorithm,  we  may  consider  a  slightly  different  network  model  that  allows  the 
state-transition  function  to  be  executed  in  multiple  rounds.  When  this  happens,  the 
message  is  generated  by  using  the  logical  state  at  the  previous  round. 

The  last  aspect  to  consider  is  the  algorithm  halting,  that  is  a  situation  such  that 
the  network  (and  therefore  each  processor)  is  in  a  idle  mode.  Such  status  can  be 
used  to  indicate  the  achievement  of  a  prescribed  task.  Formally  we  say  that  a  dis¬ 
tributed  algorithm  is  in  halting  status  if  the  logical  state  is  a  fixed  point  for  the 
state-transition  function  (that  becomes  a  self-loop)  and  no  message  (or  equivalently 
the  null  message)  is  generated  at  each  node. 
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1.4  Network  abstract  linear  programming 

In  this  section  we  define  a  network  abstract  linear  program  and  propose  novel  dis¬ 
tributed  algorithms  to  solve  it. 


1.4.1  Problem  statement 

Informally  we  can  say  that  a  network  abstract  linear  program  consists  of  three  main 
elements:  a  network,  an  abstract  linear  program  and  a  mapping  that  associates  to 
each  constraint  of  the  abstract  linear  program  a  node  of  the  network.  A  more  formal 
definition  is  the  following. 

Definition  2.  A  network  abstract  linear  program  (NALP)  is  a  tuple  (///.  (II 
consisting  of 

(i)  &  =  (/,£cmm),  a  communication  digraph; 

(ii)  (//,  fit),  an  abstract  linear  program; 

(iii)  3$  :  H  — >  I,  a  surjective  map  called  constraint  distribution  map.  □ 

The  solution  of  the  network  abstract  linear  program  is  attained  when  all  processors 
in  the  network  have  computed  a  solution  to  the  abstract  linear  program. 

Remark  2.  Our  definition  allows  for  various  versions  of  network  abstract  linear  pro¬ 
grams.  Regarding  the  constraint  distribution  map,  the  most  natural  case  to  consider 
is  when  the  constraint  distribution  map  is  bijective.  In  this  case  one  constraint  is  as¬ 
signed  to  each  node.  More  complex  distribution  laws  are  also  interesting  depending 
on  the  computation  power  and  memory  of  the  processors  in  the  network.  In  what 
follows,  we  assume  £%  to  be  bijective.  □ 


1.4.2  Distributed  algorithms 

Next  we  define  three  distributed  algorithms  that  solve  network  abstract  linear 
programs.  First,  we  describe  a  synchronous  version  that  is  well  suited  for  time- 
dependent  networks  whose  nodes  have  bounded  computation  time  and  memory,  but 
also  bounded  in-degree  or  equivalently  arbitrary  in-degree,  but  also  arbitrary  com¬ 
putation  time  and  memory.  Then  we  describe  two  variations  that  take  into  account 
the  problem  of  dealing  with  arbitrary  in-degree  versus  short  computation  time  and 
small  memory.  The  second  version  of  the  algorithm  is  suited  for  time-dependent 
networks  that  have  arbitrary  in-degree  and  bounded  computation  time,  but  are  al¬ 
lowed  to  store  arbitrarily  large  amount  of  information,  in  the  sense  that  the  number 
of  stored  messages  may  depend  on  the  number  of  nodes  of  the  network.  The  third 
algorithm  considers  the  case  of  time-independent  networks  with  arbitrary  in-degree 
and  bounded  computation  time  and  memory. 
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In  the  algorithms  we  consider  a  uniform  network  S3  with  communication  digraph 
Sf  =  (/,  £Cmm)  and  a  network  abstract  linear  program  (Sf ,  (H,  co) ,  33) .  We  assume  33 
to  be  bijective,  that  is,  the  set  of  constraints  H  has  dimension  n,  H  =  {hi,-  ■  ■  ,  hn } . 
The  combinatorial  dimension  is  8. 

Here  is  an  informal  description  of  what  we  shall  refer  to  as  the  FloodBasis  algo¬ 
rithm: 

[Informal  description]  Each  processor  has  a  logical  state  of  5  +  1  variables  taking  values 
in  H.  The  first  5  components  represent  the  current  value  of  the  basis  to  compute,  while 
the  last  element  is  the  constraint  assigned  to  that  node.  At  the  start  round  the  processor 
initializes  every  component  of  the  basis  to  its  constraint,  then,  at  each  communication  round, 
performs  the  following  tasks:  (i)  it  acquires  from  its  neighbors  (a  message  consisting  of) 
their  current  basis;  (ii)  it  executes  the  SUBEX_lp  algorithm  over  the  constraint  set  given  by 
the  collection  of  its  and  its  neighbors’  basis  and  its  constraint  (that  it  maintains  in  memory), 
thus  computing  a  new  basis;  (iii)  it  updates  its  logical  state  and  message  using  the  new  basis 
obtained  in  (ii). 

In  the  second  scenario  we  work  with  a  time-dependent  network  with  no  bounds 
on  the  in-degree  of  the  nodes  and  on  the  memory  size.  In  this  setting  the  execution  of 
the  SIJBEX  Ip  may  exceed  the  communication  round  length.  In  order  to  deal  with 
this  problem,  we  slightly  change  the  network  model  as  described  in  Section  1.3, 
so  that  each  processor  may  execute  the  state  transition  function  “asynchronously”, 
in  the  sense  that  the  time-length  of  the  execution  may  take  multiple  rounds.  If  that 
happens,  the  message  generation  function  in  each  intermediate  round  is  called  using 
the  logical  state  of  the  previous  round.  Here  is  an  informal  description  of  what  we 
shall  refer  to  as  the  FloodBcisisMulliRound  algorithm: 

[Informal  description ]  Each  processor  has  the  same  message  alphabet  and  logical  state  as  in 
FloodBasis  and  also  the  same  state  initialization.  At  each  communication  round  it  performs 
the  following  tasks:  i)  it  acquires  the  messages  from  its  in-neighbors;  ii)  if  the  execution  of 
the  SUBEX_lp  at  the  previous  round  was  over  it  starts  a  new  instance,  otherwise  it  keeps 
executing  the  one  in  progress;  iii)  if  the  execution  of  the  SUBEX_lp  ends  it  updates  the  logi¬ 
cal  state  and  runs  the  message-generation  function  with  the  new  state,  otherwise  it  generates 
the  same  message  as  in  the  previous  round. 

In  the  third  scenario  we  work  with  a  time-independent  network  with  no  bounds 
on  the  in-degree  of  the  nodes.  We  suppose  that  each  processor  has  limited  memory 
capacity,  so  that  it  can  store  at  most  D  messages.  The  memory  is  dimensioned  so  to 
guarantee  that  the  SIJBEX  lp  is  always  solvable  during  two  communication  rounds. 
The  memory  constraint  is  solved  by  processing  only  part  of  the  incoming  messages 
at  each  round  and  cycling  in  a  suitable  way  in  order  to  process  all  the  messages  in 
multiple  rounds. 

Here  is  an  informal  description  of  what  we  shall  refer  to  as  the  FloodBasisCy- 
cling  algorithm: 

[Informal  description]  The  first  S  +  1  components  of  the  logical  state  are  the  same  as  in 
FloodBasis  and  are  initialized  in  the  same  way.  A  further  component  is  added.  It  is  sim¬ 
ply  a  counter  variable  that  keeps  trace  of  the  current  round.  At  each  communication  round 
each  processor  performs  the  following  tasks:  (i)  it  acquires  from  its  neighbors  (a  message 
consisting  of)  their  current  basis;  (ii)  it  chooses  D  messages  according  to  a  scheduled  pro¬ 
tocol,  e.g.,  it  labels  its  in-neighboring  edges  with  natural  numbers  from  1  up  to  indegW 
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and  cycles  over  them  in  increasing  order;  (iii)  it  executes  the  SUBEX_lp  algorithm  over  the 
constraint  set  given  by  the  collection  of  the  D  messages  plus  its  basis  and  its  constraint  (that 
it  maintains  in  memory),  thus  computing  a  new  basis;  (iv)  it  updates  its  logical  state  and 
message  using  the  new  basis  obtained  in  (iii). 

Remark  3.  For  the  algorithm  to  converge  it  is  important  that  each  agent  keeps  in 
memory  its  constraint  and  thus  implements  the  SUBEX_lp  on  the  bases  received 
from  its  neighbors  together  with  its  constraint.  This  requirement  is  important  be¬ 
cause  of  the  following  reason:  no  element  of  a  basis  B  for  a  set  G  C  H  needs  to  be 
an  element  in  the  basis  of  GU  {/;  }  for  any  h  £  H\G.  □ 

We  are  now  ready  to  prove  the  algorithms’  correctness. 

Proposition  1  (Correctness  of  FloodBasis).  Let  .9' be  a  synchronous  time-dependent 
network  with  communication  digraph  Sf  =  (/Tcmm)  and  let  (f#  ,{H  ,a>),3S)  be  a 
network  abstract  linear  program.  If  S?  is  jointly  strongly  connected,  then  the  Flood- 
Basis  algorithm  solves  ,(H  ,(o),3§),  that  is,  in  a  finite  number  of  rounds  each 
node  acquires  a  copy  of  the  solution  of  ( H ,  Co),  i.e.,  the  basis  B  ofH. 

Proof  In  order  to  prove  correctness  of  the  algorithm,  observe,  first  of  all,  that  each 
law  at  every  node  converges  in  a  finite  number  of  steps.  In  fact,  using  axioms  from 
abstract  linear  programming  and  finiteness  of  H,  each  sequence  0)(B‘  (f )),  t  £  No, 
is  monotone  nondecreasing,  upper  bounded  and  can  assume  a  finite  number  of  val¬ 
ues.  Then  we  proceed  by  contradiction  to  prove  that  all  the  laws  converge  to  the 
same  co(B)  and  that  it  is  exactly  (o(B)  =  co{H ).  Suppose  that  for  t  >  to  >  0  all  the 
nodes  have  converged  to  their  limit  basis  and  that  there  exist  at  least  two  nodes, 
call  them  i  and  j,  such  that  ft)(Z?W(f))  =  co(B^)  co(B^)  =  co(B^(t)),  for  all 
t  >  to.  For  t  =  to  +  1,  for  every  k\  £  9f>(i),  does  not  violate  B^k,\  otherwise 
they  would  compute  a  new  basis  thus  violating  the  assumption  that  they  have  con¬ 
verged.  Using  the  same  argument  at  t  =  to  +  2,  for  every  &2  G  ^o(bi),  B^1 1  does 
not  violate  B^2-.  Notice  that  this  does  not  imply  that  Bl  does  not  violate  Bkl--,  but  it 
implies  that  a>(B^- )  <  co(B^2 1).  Iterating  this  argument  we  can  show  that  for  every 
S  >  0,  every  k  connected  to  i  in  the  graph  U^0SSf(f)  must  have  a  basis  Z?W  such 
that  co(B^)  <  co{B'[k  1).  However,  using  the  joint  connectivity  assumption,  there  ex¬ 
ists  So  >  0  such  that  LjJ0  ^°Sf  (t)  is  strongly  connected  and  therefore  i  is  connected 
to  j,  thus  showing  that  co(B^)  <  co(B^).  Repeating  the  same  argument  by  starting 
from  node  j  we  obtain  that  (o{B^  1)  <  ©(fiW),  that  implies  ®(_bM)  =  thus 

giving  the  contradiction.  Now,  the  basis  at  each  node  satisfies,  by  construction,  the 
constraints  of  that  node.  Since  the  basis  is  the  same  for  each  node,  it  satisfies  all  the 
constraints,  then  co(B)  =  co(H).  ■ 

Remark  4.  Correctness  of  the  other  two  versions  of  the  FloodBasis  algorithm  may 
be  established  along  the  same  lines.  For  example,  it  is  immediate  to  establish  that  the 
basis  at  each  node  reaches  a  constant  value  in  finite  time.  It  is  easy  to  show  that  this 
constant  value  is  the  solution  of  the  abstract  linear  program  for  the  FloodBasisMulti- 
Round  algorithm.  For  the  FloodBasisCycling  algorithm  we  note  that  the  procedure 
used  to  process  the  incoming  data  is  equivalent  to  considering  a  time-dependent 
graph  whose  edges  change  with  that  law.  □ 
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Proposition  2  (Halting  condition).  Consider  a  network  SC  with  time-independent, 
strongly  connected  digraph  (£  where  the  FloodBasis  algorithm  is  running.  Each 
processor  can  halt  the  algorithm  execution  if  the  value  of  its  basis  has  not  changed 
after  2diam((#)  +  1  communication  rounds. 

Proof.  First,  notice  that,  for  all  t  £  No  and  for  every  (i.j)  £  Ecmm, 

co(B^(t))<co(B^{t  + 1)).  (1.1) 

This  holds  by  simply  noting  that  B^(t  +  1)  is  not  violated  by  Z?M(f)  by  construc¬ 
tion  of  the  FloodBasis  algorithm.  Assume  that  node  i  satisfies  fiM(f)  =  B  for  all 
t  £  {to,- ■  ■  ,fo  +  2diam($f)},  and  pick  any  other  node  j.  Without  loss  of  generality  as¬ 
sume  that  to  =  0.  Because  of  equation  (1.1),  if  k\  £  ^4f>(i),  then  (^(B^'^l))  >  co(B ) 
and,  recursively,  if  &2  £  df>(k i),  then  co(B^2 1(2))  >  ©(/^(l))  >  Co(B).  Iterat¬ 
ing  this  argument  dist(i,y)  times,  the  node  j  satisfies  co(B^(dist(i,j)))  >  co(B). 
Now,  consider  the  out-neighbors  of  node  j.  For  every  kj  £  ./If; (j),  it  must  hold 
that  0)(B^3l(dist(i,7)  +  1))  >  Iterating  this  argument  dist(jf)  times,  the 

node  i  satisfies  ©(.BM(dist(!,y) +dist(j,i)))  >  co (B^(dist (i,j))).  In  summary,  be¬ 
cause  dist(i,  j)  +  dist(j,  i )  <  2diam(Sf),  we  know  that  (dist(i,y)  -F dist( j,  /) )  =  B 
and,  in  turn,  that 

co(B)  >  co(B^(dist{i,j)))  >  co(B). 

This  shows  that,  if  basis  i  does  not  change  for  a  duration  2diam((#)  +  1,  then  it  will 
never  change  afterwards  because  all  bases  B^\  for  j  £  {],...  ,n},  have  cost  equal 
to  co{B)  at  least  as  early  as  time  equal  to  diam(f^)  +  1.  Therefore,  node  i  can  safely 
stop  after  a  2diam(^)  +  1  duration.  ■ 
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1.5  Distributed  computation  of  the  intersection  of  convex 
polytopes  for  target  localization 

In  this  section  we  discuss  an  application  of  network  abstract  linear  programming  to 
sensor  networks,  namely  a  distributed  solution  for  target  localization.  We  consider 
a  set  of  (fixed)  sensors  {1, . . .  ,n}  deployed  on  a  plane.  These  sensors  have  to  detect 
a  target  located  at  position  x  £  ffi2.  Each  sensor  i  detects  a  region  of  the  plane, 
m\'  1  (x)  C  K2,  containing  the  target;  we  assume  that  this  region,  possibly  unbounded, 
can  be  written  as  the  intersection  of  a  finite  number  of  half-planes.  For  i£  { 1, . . . 
let  Ci  denote  the  number  of  half-planes  defining  the  sensing  region  of  sensor  i.  An 
example  scenario  with  q  =  2  for  i  €  { 1 is  illustrated  in  Figure  1 .2.  From  now 
on,  in  order  to  simplify  the  notation,  we  assume  that  c,-  =  c  for  all  i  €  { 1 , . . . , n},  so 
that  the  number  of  half-planes  (and  thus  the  number  of  constraints)  is  nc. 


Fig.  1.2  Target  localization:  set  measurements 


The  intersection  of  the  regions  detected  by  each  sensor  provides  the  best  estimate 
of  the  target  location,  M(x)  =  It  easy  to  see  that  M(x)  is  a  non 

empty  convex  set,  since  it  is  the  finite  intersection  of  convex  sets  all  containing  the 
position  x  of  the  target. 

Here,  we  are  interested  in  approximating  the  intersection  of  convex  polytopes  by 
means  of  a  “small”  number  of  halfplanes.  We  consider  the  following  approximation 
problem.  Given  a  finite  collection  of  convex  polytopes  with  nonempty  intersection, 
find  the  smallest  axis-aligned  rectangle  that  contains  the  intersection.  We  refer  this 
rectangle  as  the  “bounding  rectangle.” 

The  bounding  rectangle  has  two  important  features.  First,  the  rectangle  provides 
bounds  for  the  coordinates  of  the  target.  That  is,  let  p0  =  (pla,pl)  be  the  center  of  the 
rectangle  with  sides  of  length  a  and  b  respectively.  For  any  p  =  (p] ,  p2)  £  M(x),  then 
Ip1  —pj,\  <a/2  and  |/r  —  p2  \  <  b/2.  Second,  the  bounding  rectangle  is  characterized 
by  at  most  four  points  of  the  polytope  or,  equivalently,  by  at  most  eight  halfplanes. 
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It  can  be  easily  shown  that  computing  the  bounding  rectangle  is  equivalent  to 
solving  four  linear  programs  respectively  in  the  positive  and  negative  directions  of 
each  reference  axis.  More  formally,  let  vq  G  S 1  be  a  vector  forming  an  angle  9  with 
the  first  reference  axis.  Given  a  set  of  half-planes  H  =  {h\ .  hn},  hi  C  R2  for 
i  G  {1 , . .  • ,  n},  we  denote  (II .  (Qq  )  the  linear  program 

min  VqX 

subj.  to  afx<bi,  ie{l,...,n} 

where  hi  =  U  G  R2  |  afx  <  bt,ai  G  R2  and  bj  G  R}.  The  bounding  rectangle  may 
be  computed  by  solving  the  linear  programs  (H,C0g),  9  G  {().  n /2.  n.  3 7T / 2 } .  An 
example  is  depicted  in  Figure  1 .3 


(H, 


Fig.  1.3  Target  localization:  bounding  rectangle 


Remark  5.  (i)  A  tighter  approximation  of  the  intersection  may  be  obtained  by  choos¬ 
ing  a  finer  grid  for  the  angle  9.  Choosing  angles  at  distance  2n/k,  k  >  4,  the 
intersection  is  approximated  by  a  fc-polytope. 

(ii)  An  inner  approximation  to  a  polytope  in  dimension  d  with  np  facets  can  also 
be  computed  via  the  largest  ball  contained  in  the  polytope.  This  center  and  ra¬ 
dius  of  this  ball  are  referred  to  as  the  incenter  and  inradius  of  the  polytope.  It 
is  known  [14]  that  the  incenter  and  the  inradius  may  be  computed  by  solving  a 
linear  program  of  dimension  d+  1  with  np  constraints. 

(iii)  The  computation  of  a  bounding  k-polytope  is  a  sensor  selection  problem  in  which 

one  wants  to  select  a  few  representative  among  a  large  set  of  sensors  by  optimiz¬ 
ing  some  appropriate  criterion.  Indeed,  the  2k  sensors  solving  the  problem  are 
the  only  ones  needed  to  localize  the  target.  □ 

In  the  following  we  want  to  design  a  distributed  algorithm  running  on  a  network 
to  approximate  the  intersection  of  planar  convex  polytopes.  We  assume  the  sensor 
network  may  be  described  by  the  mathematical  model  introduced  in  Section  1.3.  Let 
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Sf  =  (/.  £wsn)  be  the  associated  communication  graph,  where  the  set  I  =  {1, . . .  ,n} 
is  the  set  of  identifiers  of  the  sensors  and  £w sn  is  the  communication  edge  map. 
We  assume  that  ^  is  a  fixed  undirected  connected  graph.  We  consider  the  network 
(abstract)  linear  programs  defined  by  (5f ,  (//,  (Og),3§)  where  (//,  a>g)  are  the  linear 
programs  defined  above  and  3$  is  the  mapping  associating  to  each  node  the  c  con¬ 
straints  describing  its  sensing  region.  Therefore,  the  distributed  algorithm  to  com¬ 
pute  the  bounding  rectangle  (k-polytope)  consists  of  4  (k)  instances  of  FloodBasis 
running  in  parallel  (one  for  each  linear  program).  We  denote  FloodRect  ( FloodPoly ) 
the  algorithm  consisting  of  the  4  ( k )  instances  of  FloodBasis  running  in  parallel  and 
the  routine  to  compute  the  bounding  rectangle  (polytope). 

Next,  we  are  interested  in  estimating  the  position  of  a  moving  target.  The  pro¬ 
posed  algorithm  can  be  generalized  according  to  the  set  membership  approach,  de¬ 
scribed  for  example  in  [11].  The  idea  is  to  track  the  target  position  by  means  of 
a  prediction  and  “measurement  update”  iteration.  The  idea  may  be  summarized  as 
follows.  We  consider  the  sensor  network  described  above,  but  with  the  objective 
of  tracking  a  moving  target.  The  sensors  can  measure  the  position  of  the  target  ev¬ 
ery  T  £  N  communication  rounds,  so  that  during  the  T  rounds  they  can  perform 
a  distributed  computation  in  order  to  improve  the  estimate  of  the  target.  The  tar¬ 
get  moves  in  the  sensing  area  with  bounded  velocity  |v|  <  vmax.  That  is,  given  its 
position  x(t),  the  position  after  T  £  N  communication  rounds  may  be  bounded  by 
x(t  +  t)  £  B(x(t):v max'!1),  where  we  have  assumed  the  inter-communication  interval 
to  be  of  unit  duration. 

In  order  to  simplify  notation  we  assume  that  each  sensor  can  measure  only  one 
halfplane  containing  the  target,  i.e.,  we  set  c  —  1.  Also,  we  suppose  that  the  sensors 
may  keep  in  memory  k  past  measures  and  use  them  to  improve  the  estimate. 

We  begin  the  algorithm  description  by  discussing  the  prediction  step.  Let  (t)  = 
{x  £  M2  |  be  the  halfplane  (containing  the  target)  measured  by 

sensor  at  time  t.  Since  the  target  moves  with  bounded  velocity,  it  follows  easily 
that  at  instant  t  +  T  the  target  will  be  contained  in  the  halfplane  h  ‘  (f  +  r|f)  =  {x  £ 
R2  |  v  maxTa^  (t)/\\a-^  (r)  1 1 )  <  b  ‘  (f)}.  The  idea  is  illustrated  in  Figure  1.4. 


Fig.  1.4  Constraint  after  the  prediction  step 


We  are  now  ready  to  describe  the  estimation  procedure.  Informally,  between 
two  measurement  instants  each  sensor  runs  a  FloodRect  algorithm  such  that  each 
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FloodBasis  instance  has  constraints  given  by  the  current  measured  halfplane  and 
the  prediction  at  time  f  of  the  latest  k  measures.  More  formally,  each  node  solves  the 
network  linear  programs  (f£ ,  (H (t) ,  (Og) , £$) ,  6  £  {0,7r/2,7r,37r/2},  where //(f)  = 
{//['] (f //[n](r)},  with 

Then  each  node  computes  the  corresponding  bounding  rectangle.  The  following 
result  follows  directly. 

Proposition  3.  Let  t  £  FI  be  a  time  instant  when  a  new  measure  arrives.  We  denote 
RecM  (; t  +  t)  the  estimate  of  the  bounding  rectangle  at  instant  t  +  T  obtained  by  run¬ 
ning  the  FloodRect  algorithm  for  the  network  linear  programs  (fS ,  (//(f),  (Oq),P§), 
0  £  {0,7r/2,7r,37r/2}.  Then 

( i )  For  any  T  £  [0,  T]  and  any  i  £  {1 , . . . , n},  the  rectangle  RecM  ( f  +  t)  is  a  subset 
of  the  rectangle  RecMft)  and  it  contains  the  target. 

( ii )  For  sufficiently  large  T,  there  exists  To  £  [0,  T]  such  that  RecM  (f  +  t)  =  RecM  (f  + 
To)  =  Recto  for  all  To  <  T  <  T,  where  Recto  solves  the  network  linear  programs. 

Proof.  To  prove  statement  i)  first  note  that  for  any  i  £  {1, . . .  ,n}  each  constraints  in 
//['l  (f )  contains  the  target.  Therefore  each  estimate  of  the  bounding  rectangle  will 
contain  the  target.  The  monotonicity  property  of  RectW  (f  +  t)  follows  by  the  mono¬ 
tonicity  property  of  each  FloodBasis  algorithm  solving  the  corresponding  network 
linear  program. 

Statement  ii)  follows  easily  by  the  fact  that  each  FloodBasis  algorithm  running 
in  parallel  solves  the  respective  network  linear  program.  ■ 


1.6  Conclusions 


In  this  paper  we  have  shown  how  to  solve  a  class  of  optimization  problems,  namely 
abstract  linear  programs,  over  a  network  in  a  distributed  way.  We  have  proposed 
distributed  algorithms  to  solve  such  problems  in  network  with  various  connectivity 
and/or  memory  constraints.  The  proposed  methodology  has  been  used  to  compute 
an  outer  approximation  of  the  intersection  of  convex  polytopes  in  a  distributed  way. 
In  particular,  we  have  shown  that  a  set  approximation  problem  of  this  sort  may  be 
posed  to  perform  cooperative  target  localization  in  sensor  networks. 
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